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We define a new unstable state in the Friedrichs model of a two-level atom. This unstable state is 
a complex eigenstate of the time evolution operator exp(-iHt) with a restricted test function space, 
which is obtained from causality conditions. The unstable state shows exact exponential decay for 
t > 0. Its emitted field is confined inside the future light-cone. In this way the long-standing problem 
of exponential catastrophe is removed. This is an example of quantum mechanics outside Hilbert 
space, which consists of generalized eigenstates in a distribution space, and a dual (test function) 
space. 
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I. INTRODUCTION 

The problem of defining unstable states has a long 
and controversial history. Although unstable particles 
abound in nature, usual quantum states with real ener- 
gies in Hilbert space describe only stable particles. Thus 
the question arises whether we can construct an unstable 
state that represents an unstable particle. 

A number of people have studied this problem. Gamow 
first introduced complex energies to model unstable par- 
ticles with exponential decay pj. Nakanishi 0] intro- 
duced complex distributions to define a complex eigen- 
state of the Hamiltonian in Lee's model 3]. The real part 
of the eigenvalue gave the particle's mass, and the imag- 
inary part gave the lifetime. In this way, a state with 
complex eigenvalue represented the unstable state. Su- 
darshan, Chiu and Gorini Q constructed complex eigen- 
states using contour deformation in the complex plane. 
Bohm and Gadella p| constructed complex eigenvectors 
using poles of the S matrix and Hardy class test functions 
(see also @). Prigogine and collaborators studied exten- 
sively the properties of complex spectral representations 
in the Friedrichs modelpj, and defined unstable states 
in Liouville space (see [8j and references therein). Still, 
the exponential growth of the field component outside 
the light cone (also called the exponential catastrophe) 
remained as a problem 

In this article we show another way of constructing 
an unstable state without exponential catastrophe in the 
Friedrichs model. This is done by separating the pole 
contribution using a suitable integration contour and test 
function space. The state we construct becomes a com- 
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plex eigenstate of the the time evolution operator e~ 
within a suitable test function space. 

The paper is organized as follows. In section II we ex- 
plain the Friedrichs model. In section III we review a 
previous approaches based on contour deformation 0, [|| 
(in Appendix A we review other approach based on a 
"rigged" Hilbert space with Hardy-class test functions 
[j|). We point out difficulties of these approaches in de- 
scribing unstable states. In section IV, we propose an- 
other way of taking the complex pole and show that this 
method eliminates the exponential growth. In section 
V, we conclude our result and discuss the extension of 
quantum mechanics outside the Hilbert space. 

II. MODEL 

We consider the Friedrichs model in one dimension 
[Til IT^ . This is a simplified version of the Lee model of 
unstable particle in the one-particle sector §. It is also 
a model of a two-level atom interacting with the electro- 
magnetic field in the dipole and rotating wave approx- 
imations 0| . Hereafter we focus on the atom-photon 
interpretation of the model. The Hamiltonian is given 
by 

/oo 
dk u k \k){k\ 
-co 

/CO 
dkv k (\l)(k\ + \k)(l\) (1) 
-oo 

where we put c = h = 1. The state |1) represents the 
bare atom in its excited level with no field present, while 
the state \k) represents a bare field mode ("photon") of 
momentum k together with the atom in its ground state. 
The excited state is analogous to an unstable particle 
state, while the photon is analogous to decay products. 
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The energy of the ground state is chosen to be zero; 
u>i is the bare energy of the excited level and Lo k = \k\ is 
the photon energy. A is a small dimensionless coupling 
constant (A <C 1). We shall consider a specific form of 
the interaction potential 



1/2 



i + (uj k /My 



(2) 



The constant M~ l determines the range of the interac- 
tion and gives an ultraviolet cutoff. Other forms of po- 
tential (form factors) may be treated in a similar way, the 
only condition being that they are exponentially bounded 
at infinity, as we will discuss later on. 

From the dispersion relation Lo k = \k\, the free- 
Hamiltonian eigenstates |fc) and | — k) have the same 
eigenvalue Lo k . We remove this degeneracy by rewriting 
the Hamiltonian 



/>oo 

ffjr=Wl|l)(l|+ / dk 

Jo 



Lo k (\S k )(S k \ + \A k )(A k \) 



+ dkV2\v k (\l)(S k \ + \S k )(l\) 
Jo 



where 



\Sk) = ^=(1*) 



\A k ) S T2 (\k) 



(3) 



*))■ (4) 



From Eq. we see that the discrete eigenstate |1) only 
interacts with the symmetric field eigenstate \S k ). Also 
in this form the Hamiltonian is expressed with energy 
eigenstates, not the k eigenstates. The anti-symmetric 
field component acts like a free field and can be treated 
separately. From now on, we concentrate on only the 
discrete atom state and the symmetric field states of the 
Hamiltonian. Changing integration variable k to lo and 
rewriting 

\Lu) = \S k ), v u = y/2v Uh , (5) 
we get the atom-field interaction Hamiltonian 



H = ux\\){\\+ / duo lo\lo){lo\ 
Jo 



du\v u {\\){u\ + \u)(l\). 



(6) 



This Hamiltonian has an exact diagonalized form. 
When the atom eigenfrequency lo\ is outside the field 
spectrum ( lo\ < 0) , we call this stable case. In this case 
the Hamiltonian is diagonalized as 



fl. = w 1 |&><&| + / dLJL0\tt)$5\ 
Jo 



where \4>i) and \cpu,) are given by 



(7) 



(8) 



\4>t) = I") 



Tf^ipS) ?? ± (w) Jq uj — uj'±ie 



dJ XV ^'\ (9) 



with 



Nx= 1 



V ± (z) 



, A V,., 

auj- 



duj- 



(10) 



(11) 



We can choose + branch or — branch for the diagonalized 
solution. These branches correspond to outgoing and in- 
coming waves, respectively. In Eq. (| 1 1 1> . l/(z ± — lo) means 
that z is analytically continued from above (+) or below 
(— ). For real z, it can be understood as 



1 



1 



lim 

e-*0+ Z ± it 



(12) 



For the stable case (lo\ < 0), the diagonalized Hamilto- 
nian has a renormalized atom state \<pi) with renormal- 
ized atom frequency lo\ < satisfying the relation 



UJ\ — LO\ 



(13) 



When the atom eigenvalue loi is inside the continuum, 
the situation changes. For 



LOl > 



xW 



(14) 



the equation ^(z) — does not have a real solution. 
Unlike stable case, we cannot maintain the renormalized 
atom state with real eigenvalue. 

One diagonalized solution for this case is due to 
Friedrichs [12|, and has the form 



H= I cLolo\F±)(F l 
'o 

, , Xv. 



±1 

LO I ' 



(15) 



ho) 



^(lo) 



dw'- 



Xv t , 



LO ■ 



■ lo' ±ie 



(16) 



The eigenstates satisfy the eigenvalue equation as well 
as the orthonormality and completeness relations 



H\F: 



= LO\F,l 



(17) 



(F±\F$) = 5(lo~lo>) 

/ dLo\F±)(F±\ = |1)<1|+ / <Lj\lo)(lo\. (18) 
Jo Jo 

Note that this solution contains only field modes. The 
bare unstable atom is viewed as a superposition of the 
field modes. The difficulty is that in this view the un- 
stable state has the memory of its creation. The decay 
law is not strictly exponential, and we can distinguish old 
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atoms and young atoms. According to this view, unsta- 
ble particles are also distinguished by their creation pro- 
cess. Because of these complications, we want another 
definition of unstable states describing indistinguishable 
particles which are independent of their creation process. 
This requires strict exponential decay with no memory 

& 

In figure^we show the survival probability of the state 
As seen in figure |2 the survival probability shows 
non-exponential decay around t = (Zeno effect). Figure 
13 shows the field generated by the initial condition of 
excited atom state and no field. We define the field bra 
(tp(x)\ as 



{if>(x)\ = / dk -— e lkx {k\ 



P(t) 




duo \ — cos(wa;)(w|, (19) 

IQ V u 

The generated field is a superposition of field associated 
with the Zeno effect, exponential field due to spontaneous 
emission and dressing cloud around the atom at x = 
[l5j| . Note that the field disappears rapidly outside the 
light cone, defined by \x\ — ct with c = 1. 





Time 



FIG. 2: The survival probability of excited atom |(l|e" i£rt |l)| 2 
near t — 0. In this short time period (Zeno time), the decay 
is not exponential. 
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FIG. 3: Fhe field intensity \(Tp(x)\e' iHt \l)\ 2 in space at t 
10. 



FIG. 1: The decay of survival probability |(l| e - ifft |l)| 2 . 

One way to get pure exponential decay is to construct 
eigenstates with complex eigenvalues. This has been 
done already, but previous constructions had their own 
difficulties, for example the exponential growth of the 
emitted field outside the light cone (exponential catastro- 
phe). In the following sections, we review the approach 
based on contour deformation. 



III. UNSTABLE STATE USING CONTOUR 
DEFORMATION 

In this section we review the construction of unsta- 
ble state through contour deformation. The construc- 



tion of complex eigenstate through contour deforma- 
tion was done by Sudarshan, Chiu, and Gorgini Q. 
Later by the perturbation expansion with regularization 
rules, Petrosky, Prigogine and Tasaki [7j investigated the 
Friedrichs model and showed that the system can be de- 
scribed as a sum of a discrete complex eigenstate plus 
continuum states. They showed that their decomposi- 
tion can be also derived by contour deformation. Let us 
follow their construction in the Friedrichs model. 

In the Friedrichs model, first we note that rj + (z) = 
has a complex root z\ = oj\ — 17 = uj\ + 0(\ 2 ) for 
< uo\ ~ 0(1). We can consider this as a complex eigen- 
value which coincides with the original discrete eigen- 
value u>i in the limit A — > 0. By contour deformation, we 
can separate z\ pole from l/r/ + (z) in the completeness 
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relation. 

/>oo 

1 = / du,\F+)(F+\ 
Jo 

= f du\F+)(F+\+ [ du\F+)(F+\ (20) 
Jr Jc 

where the contours T and C are shown in Fig0| 

A 

lm[w] 

Re[w] 

> 

r 

FIG. 4: the contours F and C 

The pole part of the contour J c can be written as 

/ ^|F+)<F+| = |<M^i| (21) 
Jc 

where 

\^=N^(\i )+ r^^m, (22) 

V Jo z{-ujJ 

(k\=N^((l\+ f^^MV (23) 

V Jo z{-uij 

{ f°° A 2 ?; 2 \ _1 

~> s H ^w^r) ■ (24 » 

This complex eigenstate |0i) satisfies the eigenvalue 
equation 

H\ih) = (25) 

Note that |0i) cannot be in the Hilbert space since it has 
a complex eigenvalue. 

It should be noted that when we act a function on \4>\) 
the function should not blow up on the deformed con- 
tour. We need suitable test functions depending on our 
choice of contour deformation. Without the test function 
consideration, unphysical growth (exponential catastro- 
phe) appears. This exponential catastrophe is a main 
difficulty of accepting this complex eigenstate as a repre- 
sentation of unstable states. 



To see this problem let us consider the time evolution 
of |(^i ) for the atom component and field component. The 
atom component of \4>i) is given by 

(l|e^ m |0i) = Nl /2 e~ lZlt . (26) 

Eq. (25 holds for all t. For t < 0, the RHS of Eq. ^ 
grows exponentially. If we had chosen the — branch of 
the Friedrichs eigenstates as a starting point, then we 
would have a similar problem: the states would decay 
for t < and would grow exponentially for t > 0. 

This exponential growth also appears in the field com- 
ponent of \4>\). The time evolution of this component is 
becomes 

(V>(x)|e- iHt |^) = e-^MsJIfc) ~ 0(e- fal <*-H>) 
for large \t- \x\\. (27) 

The field component shows exponential growth in x. The 
\(i(>(x)\e~ lHt \<j)i)\ 2 plot in x space for a fixed time is 
shown in Fig0 



field intensity 




0.02 
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FIG. 5: \(tl)(x)\e~ im \4n)\ 2 plot at t = 10. We see the expo- 
nential catastrophe for |a;| — * oo. 



We can avoid the exponential growth by choosing suit- 
able test function spaces. One such approach is due 
to Bohm and Gadella They defined complex states 
(Gamow vectors) through the poles of the S-matrix, and 
restricted their test function space to Hardy class func- 
tions from below. In Appendix A we review their ap- 
proach and also discuss difficulties in their approach. 

In the next section we propose a new way to construct 
the unstable state for the Friedrichs model. We separate 
the pole according to the type of test function, and dis- 
cuss the advantages of our construction over the previous 
constructions of unstable states. 
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IV. A NEW UNSTABLE STATE IN THE 
FRIEDRICHS MODEL 

In this section we propose a new way of defining the un- 
stable state. We want our unstable state to be memory- 
less and have no unphysical growth. To this end, we con- 
struct a complex eigenstate of the time evolution opera- 
tor e~ lHt which gives exponential decay, and a suitable 
test function space which removes unphysical growth. 



A. Complex pole and integration contour 

The complex eigenstate is related to the complex pole 
of Green's function [77+ (z)] -1 (or the pole of S-matrix 
r)~ (lu) 1 7/ + (lu)) that can be calculated by perturbation 
from the original unperturbed eigenstate. Consider the 
emission of the field by the excited atom. We focus on 
the overlap (f\e~ tHt \l) between the emitted field and a 
wave packet (/|. We restrict our attention to the case 
in which the wave packet (f\x) is square integrable and 
localized, with compact support in space representation. 

From the completeness relation of eigenstates of H , we 
have 



(/I 



iHt I 



' 11} = (f\e~ im I du>\F+){F+\l) 




(28) 



dure 



(/|FJ>(F+|1)(29) 



In Eq. (J2HI, l/rj + (u) has the pole Z\ = u>i — ij = u>\ + 
0(X 2 ) in the lower half plane. This is the pole we want 
to extract. 

Our task is to take the residue at the pole z\ in an 
integral of the form J Q dui h(u>), where h(oj) has the pole 
z\ in the lower half of the complex plane. One simple way 
to do this is making a contour which encloses the pole 
and using the residue theorem. Suppose that we make a 
counterclockwise contour C around the pole z\. If h(io) 
is analytic in C, we have 



dw h(u) 



dio 
C u 



1 



z\ 



-h(u>)(u> — Zi) = 2irihi(z) 



where hi(u>) is defined as 

hi(u>) = (lu — zy)h(u>). 



(30) 



(31) 



Note that hi(co) is analytic function inside C. 

The enclosing contour should be chosen according to 
the test function (/|. It would not be a good choice 
of contour if the test function blows up at the contour. 
Also, the test function space should be determined by 
considering underlying physics. 

In scattering experiment usually a localized wave 
packet is prepared. Say (f\x) is zero outside the region 



—xq < x < xq. In our Hamiltonian system, (f\ui) is given 

by 

(f\tj)= I dx (f\x)(x\uu) 



X ° dx(f\x) C -^l m 

-x 



According to a theorem due to Paley and Wiener [l7| . 
the function (f\co) is entire function of exponential type 
xq and belongs to L 2 on the real axis of lu (see Appendix 
iBl . This theorem shows that even though (f\ui) is L 2 on 
the positive real axis, it can be extended to the whole 
real axis and remain in L 2 . So, we can use the whole real 
axis as a part of enclosing contour. We write 



(33) 



The last term is a "background" integral, which does not 
give any pole contribution. 

To enclose the pole in the lower half plane, we need 
another piece of contour besides the real axis. If the 
function vanishes at the lower infinite semicircle, then 
the integral over the real axis is the same as the integral 
over the closed contour consisting of the real axis and the 
infinite lower semicircle, which encloses the z\ pole. For 
the Cauchy integral 



-1 
2?ri 



dco 



hx(w) 



with h\(iu) vanishing on the lower infinite semicircle, we 
can separate the z\ pole residue by subtracting the other 
pole residues. 



-1 
2^1 



dco 



E Res t 



h\(io) ■ 

jj — Z\ ' 



hi(z!) (34) 



In Eq. H34fl . p n are possible poles of h(Lu) other than Z\ 
in the lower half plane. 

The physical meaning of the test function vanishing 
at the lower infinite semicircle is causality, as we discuss 
next. 



B. Hardy-class test functions and causality 

In Eq. I)28|) we want to separate the part of the inte- 
grand that vanishes at the lower infinite semicircle in lu 
plane. This can be done through the decomposition into 
Hardy-class functions from below and from above, which 
we define now. 

A complex function G(E) on the real line is a Hardy 
class function from above (below H^) if 

(1) G(E) is the boundary value of a function G(lu) of 
complex variable (complex energy) lu = E + ir\ that is 
analytic in the half plane 77 > (77 < 0) . 

(2) 



\G{E + irf)\ 2 dE < finite 
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for all r\ with < r\ < oo (—00 < 77 < 0). 

The function G(E) is called an H± class function. 
There is an interesting relation between Hardy class func- 
tions and L 2 functions. 

Any L 2 function can be uniquely expressed as the sum 
of a function in and a function in H\_ . If a function 
f(uj) is in L 2 on the real line, we can write 



with 



We can also write 



/(w) = /+(«) + /- H 



(35) 



1 

2tt« 
-1 
2~7ri 



<2u/ 



/ + M = 



w — w — 
w' — a; + ie 



* /(t)e* 



dt f(t)e l 



&H 2 



e Hi 



(36) 



where 



/(*) 



2?T J_ 



dojf{u>)e 



f has the following properties. 



1 

-1 
2?ri 



dw 



UJ — z 



f + (z) for 7m(z) > 

for Im(z) < 

for Im(z) > 

f~(z) for im(z) < 



The Fourier transform of f ± (uj) has the property 
/+(*) = fori<0, f~(t) = fovt>0 



(37) 
(38) 

(39) 
(40) 

(41) 



The most important property for us is that a Hardy- 
class function from above (below) vanishes at the upper 
(lower) infinite semicircle in uj plane. 

From Eq. we get 



(f\e 



-iHt\ 



duj e 



1} 



77 (u) 



(/ll> 



V M Jo 



du> 



,Atv</|w'> 



uj — uj' — le 



Xv„ 



(42) 



(43) 



For t > 0, the second and third terms in the bracket [ | 
always vanish at the lower infinite semicircle in uj plane, 
if they are finite at the real line. The first term vanishes 
at the lower infinite semicircle if 



[e- iut {f\")\ 







(44) 



To see the physical meaning of this condition we write 
(see Eq. 



e-™{f\w) 



dx{f{x)e -^ c ^l (45 ) 



1 r°° 

2 V n J-oo 



Hence 



[e- iut (f\")] 
' / " dx({f\x + t) + {f\ -x-t))e™* 



(46) 



This vanishes when t > |xo|, the time required for the 
emitted outgoing wave to have an overlap with the wave 
packet. This is causality condition follows from the re- 
quirement that the integrand in Eq. 142J1 vanishes at the 
lower infinite semicircle in order to take the residue at 
the pole z\ . One point to note is that we considered the 
space and time together when we apply this condition of 
vanishing at the lower infinite semicircle. 



C. Residue at the pole and unstable state 

Now we take the residue at the pole Z\ in Eq. H42JI. 
When Eq. I|44|) is satisfied we get 

Res [</|e-^|l>] zi = (/ie-^I^X^Il) (47) 

= e-^tfl&X&ll) if [e-- t (f\uj)] + = 

where |</>i) is the complex eigenvector of Hamiltonian in 
Eq. I|22|) . To generalize this result, we define the space of 
test functions £h as the set of functions (f\e~ lHt \ui) with 
t > and (f\x) being in L 2 and having compact support. 
This also implies that (f\u>) is L 2 and exponential type 
by Paley and Wiener theorem. Due to the form factor 
of our model (see Eq. (J^J), {f\F^) is also exponen- 
tially bounded. For these test-functions we introduce a 
decomposition into a component which vanishes at the 
lower infinite semicircle in uj plane and a non-vanishing 
component. 



/H = />) + /"» 



(48) 



Precisely speaking, we denote f nv (uj) as the part whose 
maximum modulus grows exponentially as the function 
approaches the lower infinite semicircle. 

Next, we define our custom complex delta function 
8 a (uj - z\) as 



du f{u)6 a {u - z{) (49) 

duj (/» + f»)<U^ -Zt) = f v ( Zl ). 



This delta function is similar to the complex delta func- 
tions defined in Q, except that it takes only the part of 
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test functions which vanishes at the lower infinite semi- 
circle. 

If we do not restrict the test functions, we get the com- 
plex spectral decomposition 



<1lo\F+)(F+\ = I^Mil + / du;\F»d)(F+\ (50) 
o Jo 



This gives 



where 



lF +) = | W> + _^h) + r dLU > Xv ^ 



•' ' ' ' r]+ d (uj)^' ' r/+ d (uj) J co-u'+ie 

(51) 



with 



1 Z\ — LU 



■q+ d (uj) rj+(u))z+-ui 



(52) 



With delta function 8 a (u> — Z\) which restricts the test 
functions to the part which vanishes at the lower infinite 
semicircle, we devise another expression. Generally it is 
not easy to define f(co) v and f(uj) nv for function f(u>). 
But for the function (f\uj) € £n, we can define 



(f\w) v = {f\w)-, (f\u)™ = {f\u) + , 
(f\F+) v 



(53) 
(54) 



Aw w f 00 , / Afv(/|u/) 



77+ (a>) 7 w — w' + ie 



tj+(lu) J u-oj'-ie' 
-2iri\ 2 v 2 



(55) 



77+ (w) 

(/ICT (56) 

77 (w) 77 (w) Jo cj - w' - ze 

(/|F~r = (f\w)+ (57) 

Because of equations (J5^)- (|57|l , the definitions of f v and 
/"■" are different from Hardy class functions. 

With these definitions we now define the unstable state 
and its dual as 

I^WlD+f (58) 



la \=N} /2 [{l\+ I dcu^ 



- i 

00 



Zi — LU 



Zt — LU 



(59) 



where 



Zf — LU Z\ — LU 



— 2iri5 a {tu — z\). 



(60) 



(/i«^ a >=< 2 (aii) 



l la \g)=Nl /2 ((l\ 9 ) 



\v»(w\9) v 



■^ivMW , (61) 

Z\ — LU 1 



' du, 1.(62) 

Z\ — LU 



Note that if (f\Lu) nv = then (/|0 la ) = (f\fa). 

This unstable state \4>i a ) becomes a complex eigenstate 
of the e~ lHt for special kind of test functions. Using (see 
Appendix 



(K\<h*) = -^Nl /2 \ VuJ S a (Lu - zi) (63) 



and 



V -(z 1 ) = -2irt\ 2 v 2 Zi , (64) 

we have a complex eigenvalue equation (for t > 0) 

</|e-^|0 la ) =e-^ t (/|0 lQ ) (65) 

if (/|e _lfft |w) is in f-^ and vanishes at the lower infinite 
semicircle, i.e., if (/| has compact support in space and 
(e-™ t (/|a > )) + = 0. 
When (e~ tuJt (f\uu)) + ^ 0, we have 



(f\e- tHt \4>u) 

= Nl /2 e- i ^ t ((f\l)+ r cLlu 



Av w (/k) 



% - w 



+ (-2^ 1 )( e - t [e-^(/|c)+]- - (./»+)). 

(66) 

As we see later in next section the exponentially growing 
part is removed in Eq. H66JI . 

Eq. (|65|l is clearly different from the usual eigenvalue 
equation. We make comments about this equation here. 
First, this equation is an eigenvalue equation in a re- 
stricted test function space. The test function restric- 
tion is made according to the physics of the system and 
causality condition. Second, it is an eigenvalue equation 
of the time evolution operator e~ lHt , rather than H. In 
following sections we will discuss the space-time behavior 
of this new unstable state and possible complex spectral 
representations . 

Note that here we have focused on the semi-group that 
gives decay for t > 0. In a similar fashion, we can get 
results for the other semi-group with decay for t < 
starting with the — branch of the Friedrichs eigenstates, 
and exchanging the roles of functions which vanish at 
the lower infinite semicircle and functions which vanish 
at the upper infinite semicircle. 
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D. Time evolution of unstable state 

We act the time evolution operator e~ lHt on the state 
1 0i a ). The atom component of time evolved ket becomes 



(He 



-iHt 



\4>la) 



(He 



-iHt 



cLj\F-)(F-\<j> la ) 



= <1| / duj e-^lF^i^mN^XvMto - z t ) 
Jo 

= Nl /2 0(t)e- izit . (67) 

The atom component of time evolution of \<pi a ) shows 
exact exponential decay for t > 0. This is a semi-group 
time evolution. Note also the exponential growth for the 
negative t was removed. 

Similarly, the field component of \(j>i a ) is (see Appendix 

^(x)\e- iHt \cj> la ) = mx)\e-* Ht / du\F-){F~\ct> la ) 

Jo 

M l + z 2 /M 2 ' 



-2niXN} /2 ^ 77 ^ ,, „ e' tZl( ^ |a:|) 9ft- |g| 



+2AiV 1 1/2 A 



Zl 



(1 + z(/M 2 )M 



_ e -i*xt e -M\x\Q^ (68) 



-2AiVi 



1/2 



cos(w'x) 



(1 +w' 2 /M 2 )(zi +uj') 



The first term in Eq. (|68[1 comes from the complex pole 
at zi . This is the travelling field with complex frequency 
inside the light cone. It corresponds to the decay product. 
The second term and third term do not travel but decay 
with time. The second term is due to the non-locality 
of the interaction, caused by the ultraviolet cutoff in Eq. 

The third term describes the cloud surrounding the 
atom 0| . It is due to the background integral [T^| . 

None of the terms in Eq. (|68|) has exponential blowup. 
The plot of \(i/j{x)\e- tHt \(l) la )\ 2 and \(^(x)\e- lHt \l)\ 2 in 
space is shown in Fig|H| For weak coupling the field com- 
ponent of new unstable state is very close to the field 
component of bare atom decay. The field component of 
the new unstable state shows a sharp wave front, as the 
second and third terms in Eq. (|A4[) give negligible con- 
tributions. We note that if we had included virtual tran- 
sitions in the Hamiltonian, the background contribution 
would also be strictly confined within the light cone . 



E. Complex spectral representation of exp(-iHt) 

Let us apply the complex delta function 5 a (u> — Zi) 
to the complete set of Friedrichs eigenstates of iJ.The 
effect of pole enclosing contour is obtained by multiplying 
(— 2iri)(uj — zi)S a (uj — zi) to the Friedrichs solution and 
integrating over uj. The factor — 2iri appears since the 
real axis and lower infinite semicircle clockwisely enclose 
the lower half plane pole. By this operation we get the 



field intensity 
. 08 r 




FIG. 6: 

mx)\e~ 



\{^{x)\e-* Ht \^ a )\ 2 
>| 2 (thin line) at t = 



plot (thick line) and 
10. We see the field com- 



ponent of new unstable state has a sharp front at |a;j = t. 



z± residue of the part of the test function that vanishes 
at the lower infinite semicircle. 

Consider the inner product (f\e~ lHt \g) where (f\oj) 
and (w\g) are in £u- We have 



(f\e- lHt \g) = 



dtu 



dto e-*" t (/|^+)( J F+| ff )(-2«)( w - Zl )S a (w - zi) 



Jo 

x (1 - (-2ni)(w - zi)5 a (u - zi)). 



(69) 



In Eq. (|69|l . the first term in the right hand side will give 
the pole separation we wanted. By definition, S a (u> — 
zi) selects the part of e~' iU)t (f\F~) (F+\g) which vanishes 
at the lower infinite semicircle. That contains the term 
e-™(f\u)(w\g). Generally 



[e-^{f\ U ){u>\g)]- ^e-^{f\L0)-{u>\gY 



(70) 



so the first term in RHS of Eq. I|69l) is not factorizable to 
the left and right complex eigenekts. But there are cases 
in which the factorization is possible, and we show two 
simple cases. 

One case is when all terms vanish at the lower infinite 
semicircle. This condition is satisfied when t > and the 
terms [e~^ *{u\g)], [e~^ and [e~^ *</|w)<w| 5 )] 
vanish at the lower infinite semicircle (their + compo- 
nents vanish). The physical meaning of these conditions 
are the following. Suppose that (f\x) has a compact 
support [—Xf,Xf] in x space. By Eq. and Paley- 
Wiener theorem, (f\u>) is an L 2 function of exponential 
type Xf. Similarly for (x\g) with a compact support 
[— Xg,x g ], (u>\g) is an L 2 function of exponential type 
x„. When the functions approach to the lower infinite 
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semicircle, (f\w) = 0(e iux f) and (oj\g) = 0(e iux <>). So 
[e~ luJt (f\uj)} + = implies that for time t > Xf all the 
field component of (/| is inside the causal region from the 



atom. Similarly, [e 



1 + = means that for t > 



all the field component of \g) is inside the causal region 
from the atom. The condition [e~ iult {f\uj) (u;\g)} + = 
means that for time t > Xf + x g a photon from g can be 
scattered through the atom and go to the field in /. This 
is the causality condition of scattering of the field from 
the region occupied by \g) to the region occupied by (/|. 
In this case we get 



OO 



+ / due— V\F+ d ){F+\g). (71) 
Jo 

Another case that we can factorize the pole part in 
Eq. is when (f\uj) = or (uj\g) = 0. In this case 
(we consider {ui\g) = 0. Similar result can be shown for 
(f\u) = 0.), we get 



(f\e" lHt \g) 



/ z{ -u 
+ (- 2Ki\r, t )(e^ t [e^{f\w) + ]--{f\w)+)){^ 

(72) 



+ / due-™<J\F ua )(F+\g) 
Jo 

like Eq. with 



, Atv|w') 



rj +a (uj) rj +a {u>) J to — ui' + ie 

(73) 



where 



Z\ — UJ 



Eq. J72J becomes Eq. (7TJ when 



(74) 



(75) 



As explained above, this is the causality condition be- 
tween the atom and the field. When Eq. (|75|l is not 
satisfied, we have Eq. (|72|) which does not have an expo- 
nentially growing part outside the causal region. 



V. CONCLUDING REMARKS 



energy infinite semicircle) for which this pole separation 
is physically meaningful, giving a causal description of 
absorption and emission of decay products. 

The complex eigenstate constructed by separating the 
pole contribution in this way showed unique properties. 
Its atom component has exact exponential decay for the 
positive time. Its field component consists of a travelling 
wave with complex frequency inside the light cone. Thus 
the exponential catastrophe problem was removed. 

This complex eigenstate is an eigenstate of time evo- 
lution operator e~ tHt . The test function restriction is 
done by considering both time and space, rather than 
considering only time independent picture. In this way 
this complex eigenstate captures essential features of un- 
stable particles in the physically meaningful region. 

In our opinion, this work is one nice example of con- 
structing quantum mechanics outside the Hilbert space. 
The Hilbert space is very useful for describing stationary 
states. The eigenvalues are real, and the conserved norm 
represents the probability. But in our Hamiltonian, the 
atom state decays into the field, and the field is absorbed 
by the atom state. When we consider only the field space, 
norm is not conserved since fields are absorbed by the 
atom states, or the atom emits fields. In this case, we 
don't need to stick to the Hilbert space formalism, and 
distributions and suitable test functions can be used. 

In this paper the specification and decomposition of 
our test function space was done for the Friedrichs model, 
with the dispersion relation u>k — \k\. Different test func- 
tion space should be used when the dispersion relation is 
different. Moreover, we limited the initial test-functions 
to functions exponentially bounded at infinity. Inclusion 
of other functions such as Gaussians requires further con- 
sideration. Also, we have limited our discussion to Dirac 
bras or kets. An extension to density operators in Liou- 
ville space involves products of distributions, which will 
be considered in future works. 
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In this article we constructed a complex eigenstate 
with a suitable test functions which does not give an 
exponential catastrophe in the Friedrichs model. We ex- 
tended the energy spectrum to the complex plane and 
separated the complex pole. To separate the pole, we 
choose a suitable contour enclosing the pole. We estab- 
lished a class of test-functions (vanishing in the lower 



APPENDIX A: GAMOW VECTOR WITH 
HARDY CLASS TEST FUNCTIONS 

In this section we review the Gamow vector formal- 
ism introduced by Bohm and Gadella, and apply their 
formalism to the Friedrichs model. We show that the 
Gamow vector obtained also has difficulties to represent 
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the decaying state. The original derivation of Gamow 
vectors presented in this section is found in Bohm's book 

r-i 

In Ref . |16| , Gamow vectors are derived by considering 
S-matrix elements for the scattering of a pure state </> m 
into a pure physical state ip out . <fr m is a controlled free 
state and determined by the preparation apparatus. ^ out 
is a free state controlled by the registration apparatus. 

When the Hamiltonian can be written as H = Ho + V, 
where H is the free Hamiltonian and V is the interaction, 
the exact states </> + (i) and tp~(t) are written as 



if t > 0, 
- lHt if t < 0, 



(j)+{t) = (f) m {t) + / dt' G+(t - t')Vcf> in (t'), (Al) 

J — OO 

/oo 
dt' G-{t-t')V^ out {t'\k2) 
-OO 

In Eq. ljA2(l . the Green's function is given by 

(A3) 

G-(t) = [- iH> " ;^' ( A4 ) 

Defining the M0ller wave operators as 

/oo 
dt' G ± (t-t')^e-^°(*'-*', (A5) 
-OO 

Eq. (jSU and Eq. O 

can be written as 
<£+(*) = ftV n (*), tp-(t) = Q-i; out (t) (A6) 
and the scattering operator S is defined as 

s = n~^n + . (A7) 

The S-matrix element for <jf n and ip out becomes 
(^ out (t),S0 m (t)) = (O-^°"*(i),O + i "(t)) 

= (</r(i),0 + W) = (V^,0 + ) (A8) 



We can calculate Eq. (| A8|) using the eigenvectors of the 
total Hamiltonian H . If we write the eigenvectors of the 
free Hamiltonian H as \E), the eigenvectors of the total 
Hamiltonian H can be obtained using the M0ller wave 
operators 

\E ± )=n ± \E), (A9) 
H \E) = E\E), H\E ± )=E\E ± ). (A10) 

The eigenkets \E + ) and \E~) are related by 

\E+) = \E-)S(E). (All) 

Using the eigenvectors of total Hamiltonian, Eq. (|A8|I can 
be written as 



(tP-,4>+)= dE (i>-\E-)S(E)(E + \4>+) (A12) 
Jo 



We assume the S-matrix has a single complex pole Zr, 

S-l 



S(E) 



E - Z 



R 



,s + 8\(E - Z R ) 



(A13) 



To introduce the Gamow vector associated with this pole, 
Bohm and Gadella defined a test function space $_ in 
which functions are Hardy class functions from below 
and, in addition, are analytic functions that vanish faster 
than any inverse polynomial at the lower infinite semi- 
circle. It is assumed that (0 | £7 ) and (E + \(f) + ) both 
belong to <f>_. 

With these properties in mind, we continue our dis- 
cussion about the Gamow vector. When (%jj~\E~) and 
(E + \(p + ) both belong to H 2 , we have 



(^-,0+) 



dE{i]r\E-)S{E){E + \(j) + ) 



+ 



^-\E-)-^-{E+\^ 
hj — Zr 



(A14) 



Here s_i is the residue of S(E) at the complex pole Zr. 
Using Eq. (|40|) . we can write 



{<<p-,<t > +)= dE(ip-\E-)S(E)(E + \4>+) 

J —oo 

+ (-2nis- 1 )(i;-\Z R -)(Z+\cfi + ). (A15) 
Thus omitting the arbitrary vector in if£ 3 



= / dE\E+)(E^ 
Jo 

\Z^){-2-Kis^){Z+\ct> + ). 



(A16) 



where the complex eigenvector \Z R ) is given by 

^--hL^E^ (A17) 

which is a functional over Hj_ only. 

If (tjj~\E~) belongs to H^_, then the time-evolved state 



iHt \ rp- 



E ) becomes also a Hardy class function from 



below if t > 0. Using Eq. (@UJ), we obtain 

(e im i>-\Zn) = e- iZ »\^-\Z-) 

for t > and every (ip~ \E~) in H 2 _ (A18) 

This Gamow vector does show exponential decay and 
semigroup time evolution for t > 0. To see if this defini- 
tion is suitable for the representation of unstable states, 
we should also check how the field component of this 
Gamow vector behaves. We apply the above Gamow 
vector formalism to the Friedrichs model, and see how 
its field component is represented in position space. 

In the Friedrichs model, the exact eigenvectors l-E^) 
are explicitly written as \F^). From the relation 



77 + (w)-ry (uj) = 2iri\ 2 v 2 J , 



(A19) 
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we can show that 

and the scattering matrix S(u>) is 

77+ (w) 



(A20) 



(A21) 



The pole of S(u) is at Zi, which satisfies i] + (zi) = 0. 
According to the above formalism, the Gamow vector in 
the Friedrichs model is 

M = -i [°° dLJ \ F « ) 77^V (A22) 

*™ % J-00 u — Zl 

Let us calculate the field component for this Gamow 
vector. From Eqs. ljT§|) and (|A22|) . 



(A23) 



aw \/ — ; cosla; x). 

V w a; — a/ — «e 



In Eq. jA23J, we separate the component which is in 
Using Eq. (fSH|) on the second Riemann sheet to separate 
this component, we get 



A 2 HJ 



((^)l^ f »*-=2^- RV -e 



-iw|a:j 



1 1 



1 1 



A{l-iuj/M) 2 8(l-iu/M) 



(A24) 



Substituting to — z\ we obtain the Gamow vector field 
component 



(ip(x)\z 1 ) = 2ni 
1 1 



»T(*i)V 2 
1 1 



4 (1 - izi/M) 2 8 (1 - izi/M) 



(A25) 



This shows exponentially decaying behavior for x. If we 
apply e lHt ip(x) to the Gamow vector, we get 



e lHt ip(x)\ Zl ] 



2ni- 



X 2 



— iz\t 



zi fl 



1 



1 



1 



?rOi)V 2 \A{l-iz 1 /M) 2 8(l-w x /M) y 
x e -" l(t+|x|) for t > (A26) 

Since we restricted the test function space to the Hardy 
class functions from below as well as analytic functions 
that vanish faster than any inverse polynomials at the 
lower infinite semicircle (the space <f>-), the field com- 
ponent of this Gamow vector only gives the tail part of 
the exponential field e~ lZltyt+ \ x \\ which does not show 
any wavefront (see figure 0. Actually, the dominant 
part of the field emitted from the decaying atom has the 



field intensity 

O.Olr 




FIG. 7: \^{x)\e~ iHt \ Zl ){ 2 plot at t = 10. It shows the tail 
part of exponential field. 



travelling wave with wavefront at the light cone, pro- 
portional to e _lZl ^ t_ ' a; ^6'(t — The emited field term 
e -izi(t-\x\)g(j. _ originates from the e" 1 ' 1 ' term of 
(i/j(x)\F~), which is outside the space <&_. 

From the above example, we see that the emitted field 
is not fully represented by the Gamow vector \zi) with 
the test function (ip~\E~) 6 We need another 

method to separate the complex pole from the Hamilto- 
nian so that the separated pole term properly represents 
the decaying atom as well as the emitted field. 



APPENDIX B: PALEY- WIENER THEOREM 

An entire function is one which is regular for all fi- 
nite complex arguments. For the regular function f(z) 
in \z\ < \R\, we denote M(r) as the maximum modulus 
of f(z) for \z\ = r < B. For entire functions we take 
R — > 00. The entire function f(z) is called of positive 
order p and of type r if 

lim sup r~ p log M (r) = r. (0 < r < 00) (Bl) 

) >OQ 

A function of order 1 and type r (r < 00) is called a 
function of exponential type. 

The theorem by Paley and Wiener states the following. 

Theorem by Paley and Wiener. An entire function 
f(z) is of exponential type xq and belongs to L 2 on the 
real axis if and only if 



where 



f(z) = I e lz *4>{x) dx, 

-Xo 



(x) G L 2 (~oo, 00). 



(B2) 



(B3) 
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Also, if 4>(x) does not vanish almost everywhere in any xq. 
neighborhood of x (or —Xq) then f(z) is order 1 and type 



We have 



APPENDIX C: DERIVATION OF EQ. 



'■>la 



x, --> di+4^T r*s Xv "f ] y 1 i/2 fn)+ r^-^io,; 

o uj-uj' + ie) \ J 



V+ (uj) t]+{uj) 



" — i] + (uj) J (uj — uj' + ie)(z ( { — uj') 



i~ w r l + ( UJ ) V + ( UJ ) z i ~ ^ I Jo uj — uj' + ie J z+ — uj 
AT i/2( Au w \v w Xv u I ( + 
\zf — V ( UJ ) rj + (uj) z\ — uj ( 

= Nt /2 \v u (—1 — ) = 7V 1 1/2 A^(-2^)J Q (c; - Zl ). (CI) 

\ — U! Z\ — UJ J 



APPENDIX D: DERIVATION OF EQ. |K>31) 



We calculate the field component of \4>i a ) 



{i>{x)\e- iHt \4> la ) = 



1 



' , cos(u/z) ,. _ 
duj = — (u) e 



iHl 







>/2 



due - ™' 



x(-27ri)iV 1 1/2 Au w (5 a (cj - zi) 



7T 

cos(wx) 



fu V M Jo 



da/ 



dw|J£-)<F-|& a ) 
Aw w / cos(a/a;) 



1/2 



= -2mN{ 
x6 a (uj - zi 



da; 



2\e- 



■ cos(wx) 



l + w 2 /M 2 



u/(w — u' — ie) 



2A cos(u/x) 



r]-(uj)Jo (1+uj 2 /M 2 )(uj -uj' -ie) 



(Dl) 



Although (ip(x)\uj) has singularity, combined with v u the test function for 6 a (uj — z\) in Eq. i|Dl)l becomes L 2 

and analytic on the real line, grows at most exponentially at complex infinity with poles due to the form factor. Inside 
the square bracket of last term in Eq. i|Dl|) , the separation of the part which vanishes at the lower infinite semicircle 
and non-vanishing part is clear due to the exponential functions. The vanishing part inside the square bracket in Eq. 
CDlJ is 



^ e -iui(t-|a;|) 

2(1 + oj 2 /M 2 ) 



4 e(t) 



AV, 



dio 



Acos(o/x) 



r)-(w) J {1 + uj 2 /M 2 )(uj -uj 1 -ie)' 



(D2) 



Applying S a (uj — z\) to the above and rearranging terms, w e get Eq. 1)68(1 . 
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